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Abstract. The theory of diffusive acceleration of ener- 
getic particles at shock fronts assumes charged particles 
undergo spatial diffusion in a uniform magnetic field. If, 
however, the magnetic field is not uniform, but has a 
stochastic or braided structure, the transport of charged 
particles across the average direction of the field is more 
complicated. Assuming quasi-linear behaviour of the field 
lines, the particles undergo sub-diffusion on short time 
scales. We derive the propagator for such motion, which 
differs from the Gaussian form relevant for diffusion, and 
apply it to a configuration with a plane shock front whose 
normal is perpendicular to the average field direction. Ex- 
pressions are given for the acceleration time as a function 
of the diffusion coefficient of the wandering magnetic field 
lines and the spatial diffusion coefficient of the charged 
particles parallel to the local field. In addition we calcu- 
late the spatial dependence of the particle density in both 
the upstream and downstream plasmas. In contrast to the 
diffusive case, the density of particles at the shock front is 
lower than it is far downstream. This is a consequence of 
the partial trapping of particles by structures in the mag- 
netic field. As a result, the spectrum of accelerated par- 
ticles is a power-law in momentum which is steeper than 
in the diffusive case. For a phase-space density / cx p~^, 
we find s ~ Sdiff [1 + l/(2pc)], where pc is the compression 
ratio of the shock front and Sdiff is the standard result of 
diffusive acceleration: Sdiff — 3pc/(Pc — 1). A strong shock 
in a monatomic ideal gas yields a spectrum of s = 4.5. 
In the case of electrons, this corresponds to a radio syn- 
chrotron spectral index of a = 0.75. 

Key words: acceleration of particles - diffusion - plas- 
mas - shock waves - (ISM:) cosmic rays - ISM: supernova 
remnants 



1. Introduction 

Diffusive particle acceleration at shock fronts has been 
advanced as an explanation for many astrophysical phe- 
nomena involving non thermal particles (for a review see 
Blandford & Eichlcr 1987 ). A key to the success of this 
theory is its simplicity: particles are assumed to diffuse in 
space upstream and downstream of a shock front, which 
is treated as a simple discontinuity in the velocity of the 
scattering centres. The well-known result is a power-law 
spectrum of accelerated particles /(p) oc p-'^Jiff with an 
index which is independent of the diffusion coefficient, be- 
ing determined solely by the compression ratio pc of the 
shock front: Sdiff = 3pc/(Pc ~ !)• 
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In a recent Letter (Duffy et al. 1995) it was pointed 
out that the transport of particles is more complicated 
in a 'braided' magnetic field, i.e., one which has a ran- 
dom component causing field lines to wander around when 
projected onto the plane normal to their average direc- 
tion. In particular, it was shown that there may exist a 
regime of sub-diffusive transport, in which the average 
square displacement of a particle from its origin grows 
with the square-root of elapsed time. This differs sharply 
from the linear growth expected in the case of ordinary 
diffusion and reflects the restraining influence or trapping 
(in a stochastic sense) of charged particles by the magnetic 
field. At shock fronts where particles must be transported 
across the direction of the mean magnetic field in order to 
be accelerated, this leads to a modification of the acceler- 
ation rate. 

In Sect, ^we present a detailed discussion of the trans- 
port of particles in a braided field, pointing out the dif- 
fusive and sub-diffusive regimes. There are two distinct 
physical effects which come into play here: 'microscopic' 
diffusion i.e., that which would pertain in a uniform mag- 
netic field as a result of fluctuations with a length scale 
approximately equal to the gyro radius of the particle 
concerned, and 'macroscopic' diffusion of the field lines 
themselves, which results from fluctuations of much longer 
length scale. The importance of braiding, or 'field line 
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wandering' as it is called in the astrophysical literature 
(Getmantsev 1964 Jokipii & Parker 1969, Jokipii 1973) 
can be quantified in terms of a single parameter A re- 
lated to the correlation lengths of the macroscopic mag- 
netic field fluctuations and to the relative strength of the 
macroscopic and microscopic fluctuations. In the case of 
plasmas in which braiding is important, a further distinc- 
tion must be made concerning time scales. For long times, 
the transport is diffusive, and the spatial diffusion coef- 
ficient depends on a combination of the microscopic and 
macroscopic diffusivities. On short time scales, the trans- 
port is sub-diffusive. In this case, we present the analytic 
form of the particle propagator, pointing out that it tends 
to confine or trap the particles in a stochastic sense and 
is thus non-Markovian. 

These ideas are applied in Sect. || to find the distri- 
bution of energetic particles accelerated in the vicinity of 
a shock front whose normal is perpendicular to the aver- 
age direction of the field. At such shocks, the confining 
property of sub-diffusion has a strong influence on the 
acceleration process, since this is closely connected with 
the ability of particles to perform multiple crossings and 
recrossings of the shock front before being advected away 
with the downstream plasma. Assuming the angular distri- 
bution of the particles is kept close to isotropy by scatter- 
ing off the microscopic fluctuations (an assumption which 
we discuss in Sect. ||), we compute the spectrum and the 
spatial dependence of the density of accelerated particles, 
using the exact form of the sub-diffusive propagator. Sub- 
diffusion manifests itself in a steeper spectral index than 
that given by the standard formula for diffusive acceler- 
ation. It also leads to a gradient of the particle density 
in the downstream plasma - a quantity which is strictly 
uniform in the diffusive case. The transition between sub- 
diffusive behaviour on short time scales and diffusive be- 
haviour on longer ones can be investigated by constructing 
a propagator which switches from one form to the other 
at a certain point in time, which we call the decorrelation 
time. In Sect. ^ results are found using such a propagator, 
and the value of the decorrelation time for which par- 
ticle acceleration is strongly affected by sub-diffusion is 
deduced. 

Finally, in Sect. ^, we briefly discuss a few of the im- 
plications and limitations of our results. In particular, we 
point out that in synchrotron sources in which fleld line 
wandering and acceleration of electrons at shock fronts 
is important, systematically steeper spectra should be ob- 
served from locations at which the shock is predominantly 
perpendicular. 

2. Cross- field transport 



Our approach to the transport process i nvolve s separat- 
ing it into two parts (Chuvilgin & Ptuskin 1993 ). The first 
of these consists of macroscopic magnetic fluctuations of 
length scale large compared to the gyro radius of the par- 



ticle concerned, which we characterise by a relative ampli- 



tude b = {\SB\) / (B). As described in Duffy et al. ( |1995D 
we assume these fluctuations lead to a quasi-linear type 
diffusion or wandering of the field lines in the plane per- 
pendicular to the direction of the average field. This is 
described by a magnetic diffusivity Dm, which, in terms 
of cartesian coordinates with the z-axis along the average 
field, is defined by: 



2s 



D 



M 



(1) 



where Ax is the change in the x co-ordinate upon trav- 
elling a distance s along the field line. If the turbulence 
is characterised by correlation lengths and along 
and across the average field, then assuming quasi-linear 
behaviour and adopting the normalisation of Kadomtsev 
& Pogutse (|1979D, we have 



(2) 



(Note that different normalisations are used by Achterberg 
& Ball ( |1994| ) and Isichenko ( [l991a[ )). 

The second component of the transport concerns 
length scales comparable to that of the gyro radius. We 
assume that any anisotropy in the distribution function re- 
sults in the rapid growth of magnetic fluctuations (Alfven 
waves) which, in their turn, scatter the particles so as to 
remove the anisotropy. This process we model in terms 
of two spatial diffusion coefficients and responsible 
for transport along and across the local field direction, re- 
spectively. In keeping with other treatments (e.g., Achter- 



berg & Ball 1994, Jokipii 1987) we parameterise these co- 



efficients in terms of the gyro-Bohm diffusion coefficient 
Kb = ^v^mc/ {ieB), for a particle of mass m, charge e 
moving at speed with Lorentz factor 7 = (1— w^/c^)"^/^ 
in a magnetic field B: 



Kb 



CKB 

(1+^ 



(3) 



Here e < 1 is the ratio of the energy density in microscopic 
fluctuations to that in the average magnetic field. In the 
picture in which the scattering is modelled by the "t" 
operator (Chuvilgin & Ptuskin 1993), e is the ratio of the 



collision rate to the gyro-frequency. 

Thus we have on the one hand microscopic fluctuations 
of relative amplitude ^/e which are responsible for the dif- 
fusion of particles along and across the local fleld and, on 
the other, macroscopic fluctuations of relative amplitude 
b which are responsible for diffusion of the field lines. We 
shall assume e ^ 1, which means that on the microscopic 
scale, particles are closely tied to field lines along which 
they diffuse. Transport across the local field direction is 
therefore severely restricted on the microscopic scale. 



J. G. Kirk et al. 



3 



The important physical process in such a picture has 
been described by Rechester and Rosenbluth (1978): be- 



cause of the exponential divergence of neighbouring field 
lines, a particle's orbit will take it out of regions of corre- 
lated magnetic field. This can happen either because the 
particle diffuses a small distance across the field by scat- 
tering off the microscopic fluctuations, or, in the case of 
very small e, because the finite size of its gyro orbit causes 
it to encounter uncorrelated field lines (Isichenko 1991a ). 
The latter case we term "ballistic" propagation. It occurs 
if the particle decorrelates from the field before having a 
chance to diffuse microscopically, which translates into the 
condition 



e Aiir„ 



log ( — 



(4) 



(Duffy et al. |1995D . 

Ballistic propagation has been considered by Achter- 



berg & Ball (1994) in connection with the radio emission 
of supernovae. The particles propagate diffusively, with an 
effective spatial diffusion coefficient across the field given 

by 



.ballistic 



vDm 



(5) 



and the standard theory of diffusive acceleration at shocks 
applies. However, in this paper we will assume the self- 
excited microscopic turbulence is sufficiently strong to en- 
able particles to scatter before they decorrelate from the 
magnetic field. In this case, one can define a dimensionless 
parameter 



A 



(6) 



such that for A < 1 the macroscopic braiding of the field is 
irrelevant, whereas it dominates cross-field transport for 
A ^ 1. If braiding is irrelevant, there is no anomalous 
transport regime, and no modification of the diffusive ac- 
celeration picture, provided the distribution can still be 
considered isotropic (cf. Achterberg & Ball 1994). 

The most interesting parameter regime is that in which 
both braiding and microscopic scattering are important, 
which occurs for the parameter range 



^logf^l >A»1 



(7) 



Two kinds of propagation are then possible: for times less 
than that needed to decorrelate from the field id, the parti- 
cles undergo sub-diffusion, which is a combination of dif- 
fusion along a fixed field line, which itself diffuses. The 
defining characteristic of sub-diffusion is that the mean 
square deviation of a particle in the x direction from its 
position at time i = is not proportional to t as in ordi- 
nary diffusion, but rather to ^/t. 

The propagator Psuh{x,t), which is the probability of 
finding a particle in the interval (x,x + dx) at time t, given 



that it was at the origin a; = at t = 0, is found by 
simply folding the two gaussian propagators appropriate 
to diffusion of a particle along the field (i.e., in s) and of 
the field in x: 



Psuh{x,t) 



1 



+ 00 2 

ds — ;== exp 



.(8) 



(Rax & White 1992, Duffy et al. 1995). It is straightfor- 



ward to confirm that the mean square value of x increases 
as i^/^ with this propagator. Equation (8) describes the 
motion of an injected particle for times shorter than a 
decorrelation time, i.e., 



t < td 



(A^ logA)^ 

2kii 



(9) 



Subsequently, the particle decorrelates and undergoes 
"compound diffusion" , which is the coUisional transport 
regime discussed by Rechester & Rosenbluth ( 1978 ) . Here 
the propagation is diffusive in character, with an effective 
diffusion coefficient given by a combination of macro- and 
microscopic effects: 



1 



A2 



log A 



The associated gaussian propagator is 



-fcomp (3^, i) — 



: exp 



(10) 



(11) 



Although we have derived the sub-diffusive propagator 
on the basis of quasi-linear macroscopic turbulence of the 
magnetic field, it is a phenomenon which probably has a 
much wider importance. This is demonstrated, for exam- 
ple, by the numerical work of Rax & White ( 1992| ), where 
sub-diffusion is found using a simple Taylor-Chirikov map- 
ping to model the macroscopic field turbulence. Of course, 
the formulae we quote for quantities such as the decorre- 
lation time (Eq. 9) are not generally applicable to other 
models of turbulence. Furthermore, they apply only in the 
absence of significant particle drifts and only if the mag- 
netic field can be considered static. Both drifts and time- 
dependent fields may in fact be responsible for decorrelat- 
ing the particle from a field line (Isichenko |l991b ) before 
this is achieved by either diffusion across the field or by the 
finite size of the gyro orbit. Nevertheless, even though it 
does not appear possible to model the macroscopic turbu- 
lence in an astrophysical source in such detail, we may nev- 
ertheless expect the qualitative picture of a sub-diffusive 
regime on short time scales followed by a diffusive one on 
longer time scales to be a generic feature. 
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3. Acceleration at quasi-perpendicular shocks 

In a uniform magnetic field, the microscopic fluctuations 
enable particles to diffuse in space. From the theory of 
diffusive acceleration at parallel shocks, it is well-known 
that the particle distribution at the shock front is almost 
isotropic, provided that the speed of the shock is small 
compared to the speed of the individual particles. Fur- 
thermore, if the magnetic field is not directed along the 
shock normal, but is oblique, then it has been shown that 
in the absence of cross-field transport the distribution is 
almost isotropic, given that the speed of the point of inter- 
section of the magnetic field and the shock front is small 
compared to the particle speed. If the magnetic field makes 
an angle (jfnp with the shock normal as measured in the 
rest frame of the upstream plasma, this can be expressed 
as cos ^up ^ Us/v, where Ms is the speed of advance of the 
shock front into the upstream plasma and v is the par- 
ticle speed (Kirk & Heavens 1989 ). For a slowly moving 
shock front, and for relativistic particles, this condition 
is fulfilled unless the magnetic field lies almost exactly in 
the plane of the shock. In a braided field, such points are 
relatively rare. A supernova remnant, for example, has a 
shock front moving at a speed of roughly 3000 km s~^, so 
that for a relativistic electron Ug/v « 1/100 and the distri- 
bution remains close to isotropy unless the field is aligned 
to within about 1° of the shock plane. In the following, 
we assume that in a braided field, the microscopic fluctu- 
ations are able to keep the distribution close to isotropy 
at all points on the shock front and so neglect the effect 
of the small patches where the magnetic field is almost 
tangential to the shock surface. 



3.1. Sub-diffusion 

Several properties of the accelerated particle distribution 
can be found solely in terms of the single particle propaga- 
tor P(x — x' ,t — t'). In particular, given a source function 
Q{x, t) of injected particles the spatial density distribution 

n{x, t) is 



n{x, t) 



+ 00 



da;' 



dt' P{x - x' ,t ~ t')Q{x' ,t') . (12) 



Consider the case where particles are injected at a con- 
stant rate Qo for i > io on a plane boundary initially at 
x = Q which moves at constant speed u in the positive 
a;-dircction; i.e. Q{x,t) = QoS{x ~ ut)H{t ^ to). The time 
asymptotic steady state distribution, which is found by 
allowing to — oo, is 



n{y) = Qo / dtP{y + ut,t) 



where y ^ x — ut is the distance from the plane bound- 
ary. Far downstream of the boundary, y — )■ — oo, we have 
n(y) Qo/u. This can be seen by observing that far away 



from the boundary, the distribution should relax to a ho- 
mogeneous state. Seen from the rest frame of the bound- 
ary, the flux of particles across any plane y = constant 
is just un{y), and this must equal the rate of injection at 
y = 0. Alternatively, one may note that for a physically 
realistic propagator it is possible at any finite time t to 
choose spatial boundaries beyond which there are a neg- 
ligible number of particles. Denoting these boundaries by 
±a:o, we can write 



dxP(x,i) = l; P{x,t) =Ofor \x\> Xo 

I 

Then, substituting x' = x/a we have 

rxo/a 



xo/a 



dx'P{ax',t) = 1 



(14) 



(15) 



or, equivalently, P{ax,t) — s- a^^5{x), as a ^ oo. This 
property, together with Eq. (13), implies n(oo) = (far 
upstream) and n(— oo) = Qo/u (far downstream). 

On the boundary, however, the particle density de- 
pends on the precise form of the propagator. Evaluating 
the integral in Eq. (13) with the sub-diffusive propagator 
Eq. (8), one finds n(0) = 2Qo/2>U, so that the ratio of the 
density far downstream to that at the boundary is 



n(— oo) 
^(0) 



(16) 



In sub-diffusion, therefore, there is a spatial gradient 
downstream, which is absent in the case of diffusion. The 
scale length associated with this gradient can be obtained 
from the sub-diffusive propagator and is k]^'^ lU"^/^ . 

As in the case of diffusive acceleration, these prop- 
erties of the suh- diffusion- advectiori^ problem allow us 
to find the steady state phase space density fo{p) [= 
no(p)/(47rp^)] of particles accelerated at a shock, provided 
we may assume that the shock itself does not affect the 
spatial transport (see Sect. ^). Given that particles suffer 
no energy losses during propagation and are accelerated 
only on crossing the shock front, conservation of particle 
number in momentum space requires the advected flux of 
particles with momentum p far downstream of the shock 
to equal the flux of particles across this momentum level 
at the shock (assuming no particles are advected in from 
far upstream). For an almost isotropic distribution, this 
flux is simply related to the momentum derivative at the 
shock (e.g., Drury |l987| , Kirk et al. |1994| ). Writing the dis- 
tribution function far downstream as f2{p), and denoting 
by Ml and U2 the upstream and downstream flow speeds 
in the shock one finds 



^ (^iui^u,)fo{p) ] +4Vti2/2(p) = 



(13) — 



giving a power law solution /o,2 oc p with 

/2 



1 



foiPc 



1 



(17) 



(18) 
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In the case of sub-difTusive propagation, /2//0 = 
n(— oo)/n(0) = 3/2, so that 



3pc 



1 



1 



1 



(19) 



For a strong shock, pc = 4, Eq. (19) gives s = 4.5, which 
is steeper than the spectrum produced by diffusive accel- 
eration. Physically, this arises because particles are tied 
to individual field lines and are more effectively advected 
away from the shock than in the diffusive case. 

This derivation may be re-cast slightly in terms of an 
average probability per cycle, Pcscj for a particle to es- 
cape once it enters the downstream region (Bell 119781) . 
In the steady state the ratio of the far downstream flux, 
n(— oo)|m2|, to the almost isotropic flux of particles of 
speed V crossing the shock, n(0)'y/4, leads to Pose — 
[n(— oo)/n(0)](4|u2|/w). Diffusion, for which the down- 
stream gradient vanishes, leads to Pose = 4|u2|/f , whereas 
sub-diffusion gives Pcsc — Q\u2\/v- Thus, particles which 
are sub-diffusing have a higher probability of being ad- 
vected away from the shock than those which are dif- 
fusing. Assuming isotropy, the mean increase of a par- 
ticle's momentum on crossing and recrossing the shock 
is (Ap) = 4|ui — U2|p/3f from which the power law in- 
dex of the phase space distribution can be calculated: 
s — 3 + pPcsc/(Ap). For sub-diffusion, Pcsc — Q\u2\/v, 
we obtain Eq. (19). 

It is straightforward to generalise the sub-diffusive re- 
sult to any propagator of the form 



(20) 



which corresponds to transport with (Ax^) (x Us- 
ing Eq. (13), the normalisation property Ax'^{x) = 1, 
and assuming $(x) is an even function of x, one finds 
n(— 00) — n(0)(2 — (3), leading to a spectral index 



3/5c 



1-/3 



(21) 



3.2. Compound diffusion 



So far we have considered only cases in which it is a good 
approximation to use a single form of the propagator to 
describe the particle transport. However, we have argued 
that sub-diffusion applies only for times less than t^, which 
is that needed for a particle to decorrelate from its orig- 
inal field line. Especially when calculating the stationary 
density profile, it is obvious that particles injected long 
ago should really be described by the appropriate diffu- 
sive propagator Eq. (11). There are several ways in which 
a propagator could be constructed with such properties, 
and one could certainly envisage using these in a numer- 
ical simulation. However, in order keep an analytically 



^ We are grateful to Prof. Luke Drury for pointing out this 
result to us 



tractable form, we have been compelled to use a crude 
approximation. This consists of artificially changing the 
behaviour of each particle at a time id after it has been 
injected. Sub-diffusion is assumed at t < id and diffusion 
for all t > t^: The particle propagator can then be written 
as 



'Psub(a;,i) for t < td 

comp i^X ~ 2 
for t > 



^^"^'^^ " ^ /_°lclx'Pco,„p(x-a;',i-td)Psub(x',td)^^^^ 



Although this propagator permits diffusion on long time 
scales, it does not allow for any subdiffusion at all at 
t > id- In reality, one would imagine that particles perform 
a sequence of subdiffusive episodes with a mean duration 
equal to td, between which they decorrelate from the mag- 
netic structures. In contrast to this, propagator (22) sim- 
ply releases all particles from the correlated field after a 
time td- We have also investigated an alternative picture in 
which particles are released from a correlated patch once 
they have diffused along the magnetic field for more than 
a "decorrelation length" . This approach is closer to the 
physical picture of propagation in a stationary magnetic 
field, as discussed by Duffy et al. ( 1995| ). It leads to results 
which differ little from those described below. 

Introducing dimensionless variables appropriate to the 
sub-diffusive case, 



,1/3 



,,4/3 



^2/3 1/3' 



p.2/3 1/3 



(23) 



enables the dimensionless sub-diffusive and diffusive prop- 
agators to be written 



V(4r)] 



Pc 



1 pOO 1 

Psub(e,r) = - _^exp[-eV(4s) 

p(e,r) = exp [-gV(4A^compr)] (24) 



V4^ 



Here Kcomp = Hcomp/ (Dm^ kY^u'^^^) is the dimensionless 
form of the compound diffusion coefficient. This is closely 
related to the dimensionless decorrelation time Td. Eval- 
uating the mean square deviation as a function of time, 
one finds 



/oc 
-00 



for T < Td 



(25) 



(26) 



^Ay/Td/TT + 2Kcomp(T - Td) for T > Td 



Since the propagator itself is continuous at t ~ Td, so 
too is (A^^). Choosing Kcomp = 1 / V'^d means that the 
first derivative is also continuous, and this is the value 
we adopt below. (It should be noted that the arguments 
leading, for example, to Eq. (10) determine Kcomp only to 
within a factor of the order of unity.) 
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Fig. 1. Tlie particle density as a function of distance from a 
moving boundary for two values of the dimensionless decorre- 
lation time t^. The dotted line represents Mi(C,Td), i.e., those 
particles injected at the boundary within the last correlation 
time, which subsequently sub-diffused. The dashed line depicts 
'^2(^,T"d), i.e., particles injected before this. After an initial 
sub-diffusive period, these performed standard diffusion. The 
total particle density is indicated by the solid line. The length 
scale is in dimensionless units (see Eq. 23) and the density is 
normalised to unity far downstream. 



The stationary density at the shock for the mixed prop- 
agator Eq. (22) is made up of two parts: 



i(C,Td)= / drP,„b(e + T,T) 
Jo 



ni(4,Td) 

/•OO /'OO 

+ / dr / de'Pcomp(e + r-e',r-rd)Psub(e',T4^7) 

" V. ' 

n2(5,Td) 

The density ni(^, Td) consists of those particles which have 
been injected less than one decorrelation time ago, and 
have consequently propagated by pure sub-diffusion. Par- 
ticles injected before this contribute to n2(^, Td), and have 
undergone a period of Markovian diffusion after decorre- 
lating from their original field line. Evaluating these quan- 
tities involves, for ni, the numerical computation of a 
double integral. On the other hand, the expression for n2 
may be reduced analytically to a single integral. Figure [l] 
shows these quantities for two values of the dimension- 
less decorrelation time. One feature apparent from this 
figure is the irregular behaviour at a distance ^ ~ Td into 



the downstream medium. This is where, on average, our 
model propagator releases particles into diffusive propaga- 
tion. The precise form of the density profile at this point 
is sensitive to the value chosen for Kcomp. Since a realistic 
propagator would not decorrelate particles all at the same 
time, and would also reimpose sub-diffusion after decor- 
relation, the feature displayed has no direct physical sig- 
nificance. Nevertheless, this model propagator enables one 
to get a qualitative picture of the density profile further 
upstream and downstream, and also indicates the length 
scale on which a more realistic profile would vary. 




n 



4 6 



Fig. 2. The density of particles at the shock front (i.e., mov- 
ing boundary) as a function of the decorrelation time. Shown 
are the recently injected, sub-diffusing particles ni, the older, 
diffusing particles n2 and the total density. As in Fig. 0, the 
density is normalised to take the value unity far downstream. 



The power-law index of accelerated particles can be 
found from the values of ni and n2 at the shock front using 
Eq. (18). In this case, ni also reduces to a single numeri- 
cal integration. These quantities are plotted as a function 
of Td in Fig. |2| The corresponding power law index of ac- 
celerated particles is shown in Fig. ^. Here it can be seen 
that the spectrum steepens rapidly as the decorrelation 
time increases. The index associated with sub-diffusion is 
achieved for decorrelation times greater than a few in di- 
mensionless units, i.e., 



r)2/3 1/3 

,4/3 



(28) 



J. G. Kirk et al. 



7 



This quantity is approximately equal to the acceleration 



time found by Duffy et al. (1995). Thus, the steeper sub- 
diffusive spectrum applies when particles are able to in- 
crease their energy substantially before decorrelating from 
the magnetic field. 




Fig. 3. The power law index of accelerated particles at a shock 
front of compression ratio 4, computed according to Eq. (18). 
It is assumed that after injection at the shock front particles 
propagate by sub-diffusion for times less than tj. Subsequently 
they diffuse. Only if Td is short is the relatively flat spectrum 
of diffusive acceleration maintained. The time axis is in dimen- 
sionless units given in Eq. (23). 



4. Discussion 

The main conclusion of this paper is that stochastic test 
particle acceleration at a strong shock which is predom- 
inantly perpendicular does not necessarily produce the 
canonical n{p) oc spectrum expected from diffu- 

sive particle acceleration. Although much depends on the 
length scales and relative strengths of both the macro- 
scopic inhomogeneities and the microscopic scattering, 
which, taken together, determine the decorrelation time 
td, Eqs. (19) and (21) show that the spectrum produced 
depends in general on the properties of the turbulence, 
through the parameter (3 where (Aa;^) cx . 

We have made several simplifying assumptions in or- 
der to arrive at an analytic treatment of this problem. As 
in the case of diffusive acceleration, we have assumed that 
the presence of the shock front does not have an important 
effect on the spatial transport of the particles, i.e., that the 



solution of the particle transport problem at a shock front 
can be obtained from a consideration of the case in which 
an imaginary boundary moves through a uniform station- 
ary medium. This is certainly an idealisation. In a realistic 
situation the plasma upstream and downstream of a shock 
front probably supports turbulence of different character 
and amplitude. Furthermore, the plasma velocity and the 
magnetic field strength are discontinuous. Thus, we should 
allow for some change in the parameters of the propagator 
on crossing the shock, and perhaps also for a change of its 
functional form. Nevertheless, we feel that the essential 
physical feature introduced by the anomalous transport 
process - the steepening of the spectrum - is captured 
in our simple approach. This is because the nature of the 
transport in the upstream medium has very little effect on 
the particle spectrum, provided all particles entering this 
region are returned to the shock. The slope of the spec- 
trum is determined by competition between energy gain 
on crossing the shock and escape downstream. In the case 
of sub-diffusion, a particle is always tied to the same field 
line. If we consider a particle at the shock front moving 
into the downstream medium, it is evident that its escape 
probability is determined by the average distance it can 
travel along its field line before encountering the shock 
again and does not depend on the nature of that portion 
of the field line which lies upstream. The situation is dif- 
ferent if we are interested in the spatial dependence of the 
particle density upstream, or in the time taken to perform 
a cycle of crossing and recrossing, since then the nature of 
propagation in the upstream field is important. 

Another of our assumptions - that the magnetic field 
lines themselves diffuse - is not readily lifted. Our calcula- 
tion requires a specific model of the magnetic field, and we 
have too little knowledge of the properties of turbulence in 
the vicinity of shocks to do anything more realistic than 
simply assume diffusive behaviour. This corresponds to 
the case of sub-diffusion, /3 = 1/2. However, any situation 
in which the magnetic field plays a role in inhibiting par- 
ticle transport is likely to resemble that of sub-diffusion, 
or, more generally, the case (3 < 1. In this connection, it 
is interesting to note that steep spectra have been found 
in numerical simulations of acceleration in random mag- 



netic fields by Ballard & Heavens (1992). Although there 
is no well-defined average field direction in their compu- 
tations, so that the shock cannot be described as quasi- 
perpendicular, the fact that particles diffuse freely in only 
one dimension, which does not always correspond with 
the direction of the shock normal, leads one to suspect 
that here too, anomalous transport may be an important 
factor. However, as well as the effects described above, 
their results may also be infiuenced by the relativistic flow 
speeds they assumed. 

If particles are accelerated over several decades of mo- 
mentum, as is thought to be the case in supernova rem- 
nants, it may be that different types of transport dominate 
in different ranges of p. One would expect that low energy 
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particles which have relatively small gyro-radii compared 
to the correlation length of the magnetic field might be 
more closely tied to the magnetic field lines than high 
energy ones, and so suffer more from sub-diffusive type 
effects. We have not presented a detailed discussion of 
this effect here, which may, nevertheless, be important in 
the prob lem of the acceleration of cosmic rays (Duffy et 
al. |f995[ ). 

Finally, the relative steepening of the spectral index 
for /3 = 1/2 also has important implications for the accel- 
eration of electrons in astrophysical sources of synchrotron 
emission. In particular, a correlation should arise between 
the obliquity of the shock and the spectral index of the 
radiation. At perpendicular shocks, sub-diffusion can lead 
to particle spectra given by Eq. (19), which, even for un- 
cooled electrons, would produce a relatively steep syn- 
chrotron spectrum of cx i/""-^^. 
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